Scalar Glueball in Radiative J/ip Decay on Lattice 
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The form factors in the radiative decay of J/ip to a scalar glueball are studied within quenched 
lattice QCD on anisotropic lattices. The continuum extrapolation is carried out by using two 
different lattice spacings. With the results of these form factors, the partial width of J/ip radiatively 
decaying into the pure gauge scalar glueball is predicted to be 0.35(8) keV, which corresponds to 
a branching ratio of 3.8(9) x 10~ 3 . By comparing with the experiments, our results indicate that 
/o(1710) has a larger overlap with the pure gauge glueball than other related scalar mesons. 
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The existence of glueballs predicted by QCD remains 
obscure. For a scalar glueball there is some evidence of its 
existence indicated by the fact that there are ten scalar 
mesons, such as X*(1430), ao(1450) and three isoscalars 
/o(1370), /o(1500), and /o(1710). These mesons are close 
in mass and can be sorted into a SU (3) flavor nonet plus 
a glueball. Recent lattice studies predict that the lightest 
pure gauge scalar glueball has a similar mass [lH3j . Since 
the scalar glueball can mix with the nearby qq mesons, 
the three isocalars can be the different admixtures of the 
pure glueball G, the nn meson and ss meson. So the 
key problem is to identify which of the three isoscalars 
has a dominant glueball component. For this purpose, 
different mixing scenarios have been proposed by impos- 
ing different mass ordering of G, nn, and ss along with 
the known decay branching ratios of scalar mesons 0- 
[Tlj . However, the resultant mixing patterns are contro- 
versial, especially for the status assignment of /o(1500) 
and /o(1710). Obviously, more theoretical information of 
the scalar glueball is desired for the problem to be finally 
resolved. 

It is well-known that gluons can be copiously produced 
in J/ip decays because of the annihilation of the heavy 
quark pair. Among all the decays the radiative decay 
is of special importance. It is expected that the glu- 
ons produced in J/ip radiative decays dominantly form a 
glueball. If the production rate of the scalar glueball in 
the radiative decay can be reliably obtained from theo- 
retical studies, it will provide important information for 
identifying the possible candidate for the scalar glueball 
by comparing the production pattern of scalar mesons in 
these decay channels. There have been several studies on 
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this topic based on the tree-level perturbative QCD ap- 
proach and the dispersion relation method [La - tig , but 
it is difficult to estimate theoretical uncertainties in the 
used approximations. In contrast, lattice QCD provides 
the rigorous method to study the radiative decay from 
first principles. In this Letter, as an exploratory study, 
we investigate the radiative decay of J/ip into a scalar 
glueball in quenched lattice QCD. 

Gauge configurations used in this Letter are gener- 
ated using the tadpole-improved gauge action [l| on 
anisotropic lattices with the temporal lattice spacing 
much finer than the spatial one, say, £ = a s /at = 5, 
where a s and at are the spatial and temporal lattice spac- 
ing, respectively. Each configuration is separated by 500 
heat-bath updating sweeps to avoid the autocorrelation. 
The much finer lattice in the temporal direction gives a 
higher resolution to hadron correlation functions, such 
that masses of heavy particles can be tackled on rela- 
tively coarse lattices. The calculations are carried out on 
two anisotropic lattices, namely L 3 x T = 8 3 x 96 and 
12 3 x 144. The relevant input parameters are listed in Ta- 
bled! where a s values are determined from r^ 1 = 410(20) 
MeV. The spatial extension of both lattice is ~ 1.7 fm, 
whose finite volume effect was found to be small and 
negligible for glueballs Q. On the other hand, this lat- 
tice size is large enough for charmonium. For fermions 
we use the tadpole-improved clover action for anisotropic 
lattices [l7| . The parameters in the action are tuned care- 
fully by requiring that the physical dispersion relations 
of vector and pseudoscalar mesons are correctly repro- 
duced at each bare quark mass 18,]. The bare charm 
quark masses at different /? are determined by the phys- 
ical mass of J/ip, mj/^ = 3.097 GeV. The ground state 
masses of IS and IP charmonia are also calculated with 
these two lattices (see Fig. 2 and Table II of Ref. [l9| 
for the details) and the finite a s effects are found to be 
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FIG. 1: The effective energy plot for the Af + glueball with 
different spatial momenta. From top to bottom are the 
plateaus for momentum modes, p — 2nn/ L, with n — (2, 2, 2), 
(2,2,1), (2,2,0), (2,1,1), (2,1,0), (2,0,0), (1,1,1), (1,1,0), 
(1,0,0), and (0,0,0). 
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TABLE I: The input parameters for the calculation. Values 
for the coupling /3, anisotropy £, the lattice spacing a s , lattice 
size, and the number of measurements are listed. 



p 


i 


a 3 (fm) 


La s (fm) 


l 3 xT 


A^conf 


2.4 


5 


0.222(2) 


1.78 


8 6 x 96 


5000 


2.8 


5 


0.138(1) 


1.66 


12 3 x 144 


5000 



To lowest order in QED, the amplitude M for radiative 
decay J/ip — > 7G is given by 

M r , ry =e;(q,r J )(G(p f )\j»(Q)m(p i ,r)), (1) 

where q — pi — pf is the momentum of the real pho- 
ton, r and r 7 are the quantum numbers of the polar- 
izations of J/ip and the photon, respectively. e(q, r 7 ) is 
the polarization vector of the photon and j M is the elec- 
tromagnetic current operator. The hadronic matrix ele- 
ment appearing in the above equation can be obtained 
directly from lattice QCD calculation of corresponding 
three-point functions. 

One of the key issues in our calculation is to construct 
the interpolating field operator which couples dominantly 
to the so-called pure gauge scalar glueball, which is de- 
fined by using interpolating field operators built from the 
gauge fields only. For this purpose, we adopt the vari- 
ational method along with the single-link and double- 



link smearing schemes 0, Q. More specifically, since 
the irreducible representation of lattice symmetry 

group O gives the right quantum number J PC = ++ in 
the continuum limit, we construct an operator set 
{<fi a , ol = 1,2,..., 24} of 24 different gluonic operators. 
Through the Fourier transformation, 



<^a (p, t) 



E 



(2) 



we obtain the operator set {4> a (p,t), a — 1,2,..., 24} 
which couples to an + glueball state with the definite 
momentum p. For each p, by solving the generalized 
eigenvalue problem, 
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at tjj = 1, where C(t) is the correlation matrix of the 
operator set, 



(4) 



we obtain an optimal combination of the set of opera- 
tors, t) = ~Y^Va.<f>a{p}t")i which overlaps most to the 
ground state, 
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where the normalization C(p, 0) = 1 is also used. This 
is actually the case that C(t) can be well described by 
a single exponential, C(t) — We~ Et , with W usually 
deviating from one by few percents. Figure [T] shows the 
effective energy plateaus of the Af + glueball for typical 
momentum modes, where one can see that the plateaus 
start even from t = 1. 

Glueballs are noisy objects and large statistics is usu- 
ally required. In this Letter, we generated 5000 configu- 
rations for both lattice systems. In order to increase the 
statistics additionally, for each configuration we calcu- 
late T charm quark propagators Sf(x, t; 0, r) by setting 
a point source on each time slice t, which permits us to 
average over the temporal direction when calculating the 
three-point functions. Therefore, the three-point func- 
tions wc calculate in this Letter arc 
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where J M (x) = c(x)7 M c(a;) is the vector current operator, 
CVj = C7jC the conventional interpolation field for J/?/>, 
and the summation in the last equality is over all the 
possible states and vector polarizations, pi is the spatial 
momentum of the initial vector charmonium and satis- 
fies the relation pi = pf + q. The vector current J^{x), 
which is conserved in the continuum limit, is no longer 
conserved on the lattice which requires a multiplicative 
renormalization. In this Letter, we adopt the nonper- 
turbative strategy proposed by Ref. 20] to define the 
renormalization constant, 



Kit) 



4 2 l(p,t/=T/2) 



Pi 



p,t f = T/2,t) 



(7) 



(2) 

where r,,^ is the two-point function of the pseudoscalar 

(3) 

charmonium r) c , Tfjjy^ is the corresponding three point 
function with the vector current insertion on one of 
the quark lines. It should be remarked that the possi- 
ble disconnected diagrams due to the charm and quark- 
antiquark annihilation are neglected in this Letter. 

The parameters Eg, Ey, the matrix elements 
(0|$(p),0)|S(p») and {0\O v ,j\V(pi,r)) can be derived 
from the relevant two-point functions of glueballs and 
J I if). Specifically, from Eq. §5$ we have 



(0|$(p>,0)|5(p») « ^2E s (p f )V 3 . (8) 
For the vector meson we take the following convention, 
(0\OvA0)\V(p, r)> =fve j (p,r), (9) 

where fy is a parameter independent of p, and tj (p, r) 
the polarization vector of the vector meson, whose con- 
crete expression depends on reference frames and is 
irrelevant to the calculation in this Letter. By us- 
ing the multipole decomposition, the matrix elements 
(S , (p / )|J M (0)|F(p;;,r)} can be written as [H, 

E(%)l J f(°)l y ®' r )) £ j®' r ) = «w£i(Q 2 )+0mjCi(Q s 

r 

(10) 

where otnj and j3aj are known functions of pf and Pi 
(their explicit expressions are neglected here), Ei(Q 2 ) 



and Ci(Q 2 ) are the two form factors which depend only 
on Q 2 — —{pi — pf) 2 . Only the form factor E\{Q 2 ) will 
be needed to determine the decay width with 



T(J/ij> -> jG< 



o++) = t^; 



Pi 
27 M 2 



-|Si(0)| 2 , (11) 



where a is the fine structure constant, p~ i the photon mo- 
mentum with \p 7 \ = (M 2 /ip - M%)/(2Mj/$). Therefore 
we will only focus on the extraction of Ei(Q 2 ). 

To measure Ei(Q 2 ) with different Q 2 on the lattice, 
we create J / if) on lattices with the momentum pi = or 
\pi\ = 2ir/La s , and the scalar glueball with the momen- 
tum pf — 2Ttn/La s , where n is ranged from (0,0,0) to 
(2, 2, 2). Among all the combinations of the vector cur- 
rent index /i, the polarization index j, the glueball mo- 
mentum pf and the J /if) momentum pi, it is found that 
there are specific combinations which give a^i{pf,pi) = 1 
and fluiipfjPi) = 0. Hereafter, we will only select these 
combinations for our practical data analysis. An addi- 
tional benefit of this selection is that in these combina- 
tions one has Yl r e j(Pi, r ) e j(Pi, r ) = 1- 

With these prescriptions, the form factor Ei(Q 2 ) can 
be derived as, 



£i(Q 2 ,£/,i) 



4' )r(3) (p/.fl;*/.*) 

C(p f ,tf-t)T( 2 \p t ,t) 
{t,t f -t~> 0), 



ft 



2E s {p f )V 3 
(12) 



where Q 2 can be given by pi and pf, the indices of the 
three-point function L*- 3 -* and the related two-point func- 
tions F^ 2 ) are omitted here, and Zy^ is the renormal- 
ization constant of the spatial components of the vector 
current. In practice, the symmetric indices and momen- 
tum combinations which give the same Q 2 are averaged 
to increase the statistics. Traditionally, the time sepa- 
ration t and tf — t should be kept large enough for the 
ground states to contribute dominantly to the three point 
^function. Even with this large statistics, we find that 
the signal of the glueball damps rapidly with respect to 
tf —t. However, this is not a real disaster since the opti- 
mal glueball operators we use couple almost exclusively 
to the ground state, as is mentioned before. So we fix 
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FIG. 2: The extracted form factors Ei(Q 2 ) in the physical 
units. The upper panel is for /3 = 2.4 and the lower one for 
f3 — 2.8. The curves with error bands show the polynomial 
fit with Ei(Q 2 ) = Ei(0) + aQ 2 + bQ 4 , as the black dot is the 
interpolated value -Bi(O) at Q 2 — 0. 



tf — t = 1 with varying t and extract E 1 (Q 2 ) from the 
plateaus of E\ (Q 2 , t / , t ) . With the very high statistics in 
this Letter, the hadron parameters, such as the energies 
of the glueball and J/ip, the constant fy in Eq. © and 
the matrix elements (0|$(p, 0)\S(p)) can be determined 
very precisely and are treated as known parameters. 

Ei(Q 2 ) for different Q 2 are extracted from the same 
configuration ensemble and are therefore highly corre- 
lated. In the data analysis we fit them through the cor- 
related data fitting. For each lattice system, the 5000 
configurations are divided into 100 bins with 50 config- 
urations in each bin. The measurements in each bin are 
averaged and the average is taken as an independent mea- 
surement. After that, all Ei(Q 2 )s are extracted simulta- 
neously through the jackknife method. In order to get 
the form factor at Q 2 = 0, we carry out a correlated 
polynomial fit to the E^Q 2 ) from Q 2 = -1.0 GeV 2 to 
2.5 GeV 2 , 



E 1 {Q 2 ) = E 1 (0) + aQ 2 + bQ 4 



(13) 



Figure [2] shows the final results of E^Q 2 ) for j3 = 2.4 
(left panel) and (3 = 2.8 (right panel), where the red 
points are the calculated value with jackknife errors, and 
the red curves are the polynomial fit with jackknife error 
bands, the black points label the interpolated i?i(0,a s ). 

The last step is the continuum extrapolation using the 
two lattice systems. The continuum limit of £4(0, a s ) is 
determined to be -Ei(O) = 0.0536(57) GeV by performing 
a linear extrapolation in a 2 . For the continuum value of 
the scalar glueball mass, we take Mq — 1.710(90) GeV 



TABLE II: Listed in the table are the Af + glueball masses 
Mg, the renormalization constants Zy\a a ) of the spa- 
tial component of the vector current, and the form factors 
Ei(Q 2 = 0, a s ) calculated on the two lattices with /3 = 2.4 
and 2.8, respectively. Also shown are the continuum extrap- 
olation of -Ei(O) and the resultant partial width T. 



p 


M G (GeV) 


Z { y\a s ) 


£i(0,a s ) (GeV) 


T(keV) 


2.4 


1.360(9) 


1.39(2) 


0.0708(43) 




2.8 


1.537(7) 


1.11(1) 


0.0602(31) 




oo 


1.710(90) [3] 




0.0536(57) 


0.35(8) 



from Ref. Q. Thus, according to Eq. (fTTj) . we finally get 
the decay width T(J/ip -> 7G0++) = 0.35(8) keV. Using 
the reported total width of J/ip, T tot = 92.9(2.8) keV, 
the corresponding branching ratio is 



T(J/ib -> 7G 0++ )/T t ot = 3.8(9) x 10 



-3 



(14) 



By comparing our result with their production rates in 
the radiative decay of J/ip, we can get some useful infor- 
mation for the glueball components of the scalar mesons 
/ (1710), / (1500), and / (1370). From PDG2010 [Hj], 
the branching ratios of the observed radiative decay 
modes of J/ip to / (1710) are: Br{J/^ -> 7 / (1710) ^> 
^KK) = 8.5lJ-jxlO- 4 , Br(J/t[j -> 7 / (1710) -> jirn) = 
(4.0 ± 1.0) x 10~ 4 , Br(J/%l> -> 7/o(1710) -)• jujuj) = 
(3.1 ± 1.0) x 10~ 4 , which add up to about 1.5 x 10~ 3 . 
With the measured branching ratio -Br(/o(1710) — > 
KK) = 0.36 ± 0.12 jp, and the ratio r(/ (1710) -> 
7T7r)/r(/ (1710) -> KK) = 0.4lt°;^ [H (Ref. [13 also 
predicts this ratio to be 0.32 ± 0.14 from a coupled chan- 
nel study of meson-meson S- waves), one can estimate 
the production rate of / (1710) to be (2.4 ± 0.8) x 10~ 3 
or (2.7 ± 1.3) x 10 -3 . This is compatible with our lat- 
tice result. For the / (1500), PDG2010 gives a lower 
bound to its production rate in J /tp radiative decay, 
Br{J/^ -> 7/o(1500)) > 5.7(8) x lO" 4 (p. On the other 
hand, with the BESII result Br(J/ip ->■ 7/ (1500) -> 
7 7T7r) = (1.01 ± 0.32) x 10" 4 [13, and the branching ra- 
tio Br(/ (1500) -)• tttt) = 0.349 ± 0.023 [U, Br(J/i[> -)• 
7/o(1500)) is estimated to be 2.9(9) x 10" 4 . Both are 
much smaller than our prediction. Finally, there is no 
evidence of the production of / (1370) in the Jji\> radia- 
tive decays. Based on this comparison, /o(1710) seems to 
have scalar glueballs as its dominant components, while 
for the other two scalar mesons, this does not seem to be 
the case. 

To summarize, we have carried out the first lattice 
study on the E\ amplitude of J/tp radiatively decays 
into the pure gauge scalar glueball G Q ++ in the quenched 
approximation. With two different lattice spacings, the 
decay amplitude is extrapolated to the continuum limit 
with a value E X {Q 2 = 0) = 0.0536(57) GeV. Thus, the 
partial decay width T(J/ijj — > 7G0++) is predicted to 
be 0.35(8) keV, which gives the branching ratio r/r to t = 
3.8(9) x 10~ 3 . We admit that the systematic uncertainties 
due to the quenched approximation are not under control 
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in this Letter, however, we are pleased to see that a re- 
cent 2 + 1-flavor dynamical lattice study on the glueball 
spectrum claims that there are not large unquenching ef- 
fects observed, especially for the scalar and tensor glue- 
balls [24j. Anyway, our results are helpful in the sense 
that /o(1710) appears to have the largest overlap to the 
pure gauge glueball among the relevant scalar mesons. 
We hope this fact sheds some light on the long-lasting 



puzzle of the identification of the scalar glueball. 

We acknowledge fruitful discussions with Keh-Fei Liu 
and Hai-Yang Cheng. This work is supported in part 
by the National Science Foundation of China (NSFC) 
under Grants No.10835002, No.11075167, No.11021092, 
No. 10675101, No. 10947007, and No.10975076. Y. C 
and C. L also acknowledge the support of NSFC and 
DFG (CRC110). 



[1] C.J. Morningstar and M. Peardon, Phys. Rev. D 56, 4043 
(1997). 

[2] C.J. Morningstar and M. Peardon, Phys. Rev. D 60, 

034509 (1999). 
[3] Y. Chen et al, Phys. Rev. D 73, 014516 (2006). 
[4] C. Amsler and F.E. Close, Phys. Lett. B 353, 385 (1995). 
[5] C. Amsler and F.E. Close, Phys. Rev. D 53, 295 (1996). 
[6] F.E. Close and A. Kirk, Eur. Phys. J. C 21, 531 (2001). 
[7] F.E. Close and Q. Zhao, Phys. Rev. D 71, 094022 (2005). 
[8] W. Lee and D. Weingarten, Phys. Rev. D 61, 014015 

(1999). 

[9] D.-M. Li, H. Yu, Q.-X. Shen, Commun. Theor. Phys. 34, 
507 (2000). 

[10] F. Giacosa, Th. Gutsche, V.E. Lyubovitskij and A. 

Faessler, Phys. Rev. D 72, 094006 (2005). 
[11] H.Y. Cheng, C.K. Chua, and K.F. Liu, Phys. Rev. D 74, 

094005 (2006). 

[12] V.A. Novikov, M.A. Shifman, A.I. Vainshtein and V.I. 

Zakharov, Nucl. Phys. B165, 67 (1980). 
[13] B.A. Li and Q.X. Shen, Phys. Lett. 126B, 125 (1983). 
[14] J.F. Donoghue, K. Johnson and B.A. Li, Phys. Lett. 

99B, 416 (1981). 



[15] Y. Liang, K.F. Liu, B.A. Li, S.J. Dong, and K. Ishikawa, 

Phys. Lett. B 307, 375 (1993). 
[16] H.B. Meyer, J. High Energy Phys. 01 (2009), 071. 
[17] C. Liu, J. Zhang, Y. Chen, J.P. Ma, Nucl. Phys. B 624, 

360 (2002). 

[18] S. Su, L. Liu, X. Li, and C. Liu, Int. J. Mod. Phys. A 21, 

1015 (2006), Chin. Phys. Lett. 22, 2198 (2005). 
[19] Y.-B. Yang, Y. Chen, L.-C. Gui, C. Liu, Y.-B. Liu, Z. 

Liu, J.-P. Ma, and J.-B. Zhang (CLQCD Collaboration), 

Phys. Rev. D 87, 014501 (2013). 
[20] J.J. Dudek, R.G. Edwards, and D.G. Richards, Phys. 

Rev. D 73, 074507 (2006). 
[21] K. Nakamura (Particle Data Group), J. Phys. G 37, 

075021 (2010). 

[22] M. Albaladejo and J.A. Oiler, Phys. Rev. Lett. 101, 
252002 (2008). 

[23] M. Ablikim et aZ.(BES Collaboration), Phys. Lett. B 

642, 441 (2006). 
[24] E. Gregory, A. Irving, B. Lucini, C. McNeile, A. Rago, 

C. Richards, and E. Rinaldi, J. High Energy Phys. 10 

(2012), 170. 



